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Finite-temperature Drude weight within the anisotropic Heisenberg chain 
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Finite-temperature Drude weight (spin stiffness) D{T) is evaluated within the anisotropic spin-1/2 Heisen- 
berg model on a chain using the exact diagonalization for small systems. It is shown that odd-side chains allow 
for more reliable scaling and results, in particular if one takes into account corrections due to low-frequency 
finite-size anomalies. At high T and zero magnetization D is shown to scale to zero approaching the isotropic 
point A = 1. On the other hand, for A > 2 at all magnetizations D is nearly exhausted with the overlap with 
the conserved energy current. Results for the T- variation D{T) are also presented. 
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I. INTRODUCTION 



It has by now become evident that many-body (MB) quan- 
tum systems of interacting particles behave with respect to 
transport quite differently if they are either integrable or non- 
integrable 111 |2|]. In integrable systems the anomalous re- 
sponse shows up in a possibility of finite-temperature stiff- 
ness (Drude weight) D{T) > lH], both the charge (or 
spin) and the thermal one 15|], indicating the dissipationless 
d.c. transport at T > 0. The prototype model for this phe- 
nomenon is the anisotropic spin-1/2 Heisenberg model on 
a chain, equivalent to the one-dimensional (ID) t-V model 
of spinless fermions with nearest neighbor repulsion. Within 
this model one of the conserved quantities is the energy cur- 
rent JE leading to the singular but trivial thermal dynamical- 
conductivity linear response ||5|], i.e., k(w) = DtS{uj). On the 
other hand, the spin current j and the corresponding dynam- 
ical spin conductivity (diffusivity) a{uj) at T > is still the 
subject of very active theoretical investigations and debate. 

In the case of a nonvanishing projection of the spin current 
j on local conserved quantities Q„ the Mazur inequality offers 
a firm proof of finite D{T ^ 0) > |51] in the thermodynamic 
limit . Still, at zero magnetization, i.e., at the total spin S^ — 
the overlap with all Qn vanishes independent of the anisotropy 
A ||5|]. To employ the same argument one possible path is to 
construct more general nonlocal conserved quantities ||6l \M 
which should be further explored. 

The (original) alternative formulation via the MB level dy- 
namics induced in a ID ring via an external flux p, 4] offers 
a qualitative understanding and is the starting point for nu- 
merical calculations using the full exact diagonalization (ED) 
method [8-111. The latter so far did not eliminate disagree- 
ment on several questions : a) is Z? a monotonous function of 
A at fixed 5"= li^, b) does D{T > 0] vanish on approaching 
the isotropic point A = 1, S*^ = |l9lll2|l, c) which if any an- 
alytical result, obtained via the Thermodynamic Bethe Ansatz 
OISUlSll . is correct and compatible with numerical investiga- 
tions. 

In the following we present results of the numerical study 
for D{T) as obtained using the ED and the scaling for small 
systems. In contrast to previous works ||9l |lO(] we perform 



the study within the canonical ensemble which offers much 
faster convergence with the chain size L, at least approaching 
the isotropic point A ^ 1, 5^ ^ 0. To avoid quite singular 
behavior of even-lengths chains, we study spin systems with 
odd L. In particular, we pay the attention to possible low- 
frequency contributions in the dynamical conductivity a{uj) 
which can give an insight into anomalies around commensu- 
rate A = cos(7r/i') with integer i^, e.g., at A < 0.5. 

The paper is organized as follows: In Sec. II we present 
the model and Drude weight D as zero frequency contribu- 
tion to dynamical conductivity. We shortly also describe nu- 
merical method used to analyse it. Our results are presented 
in Sec. III. First we investigate the high-temperature limit 
C = TD, where we emphasize the low-frequency contribu- 
tions which can mask the correct result. We show also that 
within Ising-type regime A > 1 the Drude weight calculated 
via the overlap with the conserved energy current gives nearly 
perfect results. Finally we focus on the temperature variation 
of i:>(T). 



II. DRUDE WEIGHT 

We study the anisotropic 5 = 1/2 Heisenberg model on a 
chain with L sites and periodic boundary conditions 






(1) 



where Sf are component of the 5 = 1/2 spin operators. In 
order to define the Drude weight (spin stiffness) D it is con- 
venient to map the model ([l]i via the Jordan- Wigner transfor- 
mation onto the t-V model of interacting spinless fermions 
ad ding a fictitious magnetic flux $ = i0 through the ring 
041 Im], entering the hopping matrix elements, 

B = t5](e^*4c.+i +h.c.) + ^5] (n,- 1) (n,+i - i), 

i i 

(2) 
Ui = clci, t = J/2 and V = 2iA. Here we consider only 
chains with odd number of fermions N to avoid additional 



boundary fermionic sign and other finite-size effects discussed 
in more detail below. In the following we use everywhere J — 
1 in order to facilitate the comparison with the majority of 
previous works and references Il9l ll0lll3ll . Note that relevant 
parameters are now the total spin Sz and magnetization s = 
Sz/L or the fermion density or band filling n = N/L = s + 
1/2. 

Via the corresponding spin (particle or charge within the 
fermionic model) current 



j = f^(je'*4c,+i+h.c.), 



(3) 



one can express the dynamical (spin) conductivity at general 
temperature T > as 

a{oj) = 2ttD5{uj) + cr,.eg(w), (4) 

where the regular part Cregioj) expressed in terms of eigen- 
states \n) and eigenenergies e„, 

TT 1 - e-^" 

CTregK^) = 



L oj 



^ p„|(n|j|m)p(5(e„-e„ -w), 



e„#c„ 



(5) 

while the dissipationless component with the Drude weight 
(spin stiffness) D can be related to the flux dependence of 
MB states (|3|], in analogy with the original formulation by 
KohnlfF 



c^i&.. 



a2e„(^) 



(6) 



where pn — exp(— /3e)/Z are corresponding Boltzmann fac- 
tors. 

The relation ^ is convenient for the ED numerical eval- 
uation of D{T) in small systems, since it only requires the 
calculation of eigenvalues e„((/)). Finally we are interested in 
the result within the thermodynamic limit L — >^ oo at fixed 
T and magnetization s (filling n) , hence several strategies 
to obtain the thermodynamic value are possible. Since we 
mostly consider the high-T limit (allowing for most accurate 
ED results in small systems) and ED sizes are quite limited 
L < 21, we perform the canonical calculation at total spin S^ 
(fermion number N). The grand canonical evaluation at avail- 
able L and high T has a very broad distribution of N, leading 
to overestimates of D (or at least its slow convergence with 
L) in the vicinity of the isotropic phase, i.e., at s ^ 0, A ^ 1. 
On the other hand, also results with even L show deficiencies 
llOfl. Treating in Eq.® the flux cj) as parameter, correspond- 
ing D{(j>, T ^ 0) show strong anomaly at (/> — ^ for even L 
and even N due to the particle-hole symmetry and degeneracy 
of MB levels. In addition, even-L systems give at odd N con- 
siderably lower values for D at A < 1 and small L 11 Oil (an 
origin could be also particle-hole symmetry absent at odd L) 
remedied presumably only at much larger L. To avoid these 
complications, we in the following consider only systems with 
odd L = 5 — 21 (for L = 21 only one fc-vector due to very 
high CPU requirements) which reveal much weaker and more 
regular D((/)) dependence. 
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Figure 1. (Color online) High-T Drude weight C = TD vs. 1/L for 
zero magnetization s — and different A as obtained for systems 
with odd L = 5 - 19. 



III. RESULTS 
A. High-Temperature Limit. 

In the following we mostly concentrate on the limit 
T — > oo, expecting that obtained results are quite generic and 
qualitatively similar at any T > 0. Since for T — > oo, D{T) 
scales as l/T the relevant and nontrivial quantity is C = 
TD{T), representing also the limiting value of the current- 
current correlation function C ~ Cjj (t -^ oo) i^. Let us first 
consider the most delicate zero-magnetization s = (half- 
filling n — 1/2) case. Since we choose odd L, the actual 
calculations are performed for closest odd N — {L ± l)/2. 
Results for C vs. 1/L for all odd L ~ 5 — 19 are presented 
for different A in Fig. 1. Several conclusions can be drawn 
directly from obtained results: a) Both values as well as the 
scaling with L are qualitatively different between A > 1 and 
A < 1. It is evident that for A > 1 the only consistent limit 
appears to be C = 0. b) There are some visible anomalies 
near A < 0.5 which indicate on a nonuniform dependence of 
C(A) 1 10] and in particular different scaling L ^^ oo which 
we discuss in more detail below. 

In order to resolve the origin of the deviations of C at 
A < 0.5 as well as of quite regular convergence of results for 
other values of A we investigate the dynamical a{uj), shown 
conveniently also in the integrated form for T — > oo. 



T 

I{U)) ^ C -\ / (7reg{^')duj' 

TT 



consistent with the sum rule 



I{lo 



-^ ^kv 



-T{Hkin)/L, 



(7) 



(8) 



where Hkm is the kinetic-energy part in the model (|2]i. Ckin 
can be evaluated exactly in the /? — > limit, even for finite L 
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Figure 2. (Color online) Regular part of dynamical conductivity 



aregi'j-') and the integrated one /(a;) (inset) for s • 
A = 0.5, and b) A = 0.86 for different sizes L = 13 



and fixed N, 



^ki'. 



' A L\ L 



N -1 



1 



: and: a) 
21. 



(9) 



In Fig. 2 we present characteristic results for crreg(w) as 
well as /(cli) (in inset) for two commensurate values v — ?,,&, 
i.e., A = 0.5, a/3/2 = 0.866, respectively. We note that for 
A = 0.5 the incoherent part in CTreg(w) is quite L-independent 
in a broad range L = 13 — 21 and consequently the conver- 
gence of obtained Drude weight C vs. \/L is very stable. Less 
obvious case is A = 0.866 {v — 6) being already closer to the 
critical value A = 1. The incoherent a^egi'-^) reveals here a 
low-w contribution whereby the peak is shiftirig as with 1/L 
as observed even more pronounced for A > 1 il6ll . However, 
in the present case the peak intensity as well diminishes with 
L (a closer inspection reveals that the peak ujp also vanishes 
here faster than 1/L) so that the integrated /(w) in Fig. 2b 
appears to have well defined limit C — I{uj —> 0). 

In Fig. 3 we present /(cj) for A = 0.25 characteristic 
for the regime A < 0.5. We note that the high-w part is 
quite L-independent (note that for L=21 we calculate only 
one fc-vector, which influences slightly the sum rule /(w — > 
oo)) similar to results for A — 0.5 in Fig. 2a. However, 
there is also a well visible anomalous low-cj contribution at 
0.02 < a; < 0.08 (see the inset). The peak in o'(cj) (as 




Figure 3. (Color online) Integrated dynamical conductivity /(a;) for 
A — 0.25, s = and various sizes L. The inset focuses on the 
low-tt; regime. 



obtained from I{uj) in the inset of Fig. 3) appears to shift 
towards ujp — Q somewhat faster than 1/L (approximate fit 
ujp ~ 1.342/L — 0.017) whereas its weight in I{uj) increases 
with the system size. This deviation can be counted as an 
additional contribution to effective 6C. This is, e.g., in con- 
trast to case A = 0.866, where the intensity decreases with L 
(Fig. 2). Although the origin of the low-cj anomaly is not well 
understood it seems that it is absent for commensurate values 
of A = cos(7r/j/) which possess additional degeneracies 11 1 3r . 



Results for C vs. 1/L as in Fig. 1 can be used to extrapo- 
lated to the thermodynamic value C where we use the extrap- 
olation C(L) = C + a/L + C/L^. Obtained results for C(A) 
are presented in Fig. 4. On the other hand, one can coiTect 
C(L)with the low-w contribution (7(L) = C{L) + SC{L) and 
get modified extrapolation C, also presented in Fig. 4. We can 
now compare the results with the analytical result obtained via 
Thermodynamic Bethe Ansatz (TBA) lll3lll5ll . 



C = 



7 - ^ sin(27) 
I67 



A — 003(7), 



(10) 



the validity of which has been still questioned 0121 II5I1 . We 
note that the agreement of the analytical form ( fTOb with the 
corrected numerical C is very satisfactory for s = within 
the whole regime of A. 

Let us now turn to the dependence of C on magnetization s 
(filling n). It is evident that one gets C = within the Ising- 
type regime A > 1 only for s = 0. Results for C at A = 1.7 
and A = 3 are shown in Fig. 5 for fixed L = 19 and all avail- 
able Sz- It is indicative that C{s) are nearly equal for both 
A > 1. To go beyond the finite-size results one can also per- 
form the scaling to L — > 00 analogous to n = 1/2 case which 
is possible, e.g., for s = 1/4 (taking into account results for 
L = 5 - 25) and s == 1/3 (with results for L = 9 - 21). 
Corresponding results for the extrapolated C are also plotted 
in Fig. 5, confirming that C{s) become essentially universal 
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Figure 4. (Color online) High-T Drude weight C = TD vs. A at 
magnetization s = obtained: using the ED and finite-size scaling 
of D{L) (full curve), adding the low-o; correction (dashed curve), 
and within the analytical TBA L13; . il5,1 (dotted curve). 



Figure 5. (Color online) High-T Drude weight C vs. magnetization s 
within the Ising-type regime A — 1.7, 3. obtained for fixed L = 19, 
with the 1/L — >■ extrapolation for s = 0, 1/4, 1/3 (dots), and the 
analytical approximation C3, Eg.llllt. for A = 3. 



for A > 1. 

It has been already observed in Ref.||5|] that within the Ising 
regime A > 1 the Drude weight can be via the Mazur inequal- 
ity well exhausted with the overlap onto the simplest nontriv- 
ial local conserved quantity Q^ — je representing the energy 
current. At T -> 00 this overlap can be evaluated exactly 
leading to 



C.= 



1 (JQs)' 


A2s2(l_4s2) 


2L {QD 


l + 2A2(l+4s2) 



(11) 



From Fig. 5 we see that the agreement between the approx- 
imate C3, Eq.dTTTl. and the extrapolated C is nearly perfect 
for large A ;:^ 1, e.g., A = 3, while for A = 1.7 the value 
C3 starts to decrease, so that C3 < C. In fact we observe 
from Eq. (fTTl i that C3 just saturates as a function of A for 
A > 1.7—2 and its value there can already reasonably re- 
produce C. We stress again completely different behavior is 
for A < 1 and s = 0. In this case one gets C3 — (as well 
as higher overlaps C,i>3 = due to particle-hole symmetry), 
hence the Mazur inequality with local conserved quantities is 
unable to reproduce C > at s = |5]. 

Let us further consider the normalized Drude weight D* = 
D/ckin which represents the relative weight of the dissipa- 
tionless transport within the whole sum rule, Eq.dSJ, i.e., we 
have < D* < 1. Since one cannot perform a system- 
atic extrapolation L — > 00 for arbitrary magnetization s we 
present in Fig. 6 results for D* within the whole (half) plane 
A, s > as calculated in systems with fixed L = 19. Apart 
from some anomalies observed (without the correction 6C) 
akeady in Fig. 3 we confirm quite regular dependence D* on 
(A,s). It is quite evident that in the limiting case A = 
(XY model) we get D* — 1 corresponding to noninteract- 
ing fermions where the whole sum rule is within the Drude 
weight. The same hold for maximal magnetization s -^ ±1/2 
(for nearly empty or full band, n — > 0, n — > 1, respectively) 




Figure 6. (Color online) Normalized Drude weight D* within the 
plane A, s as calculated in systems with fixed size L = 19. 



where the interaction does not play a role. For fixed A the 
minimum of D* is always at s = whereby the dependence 
D*{s) is nearly universal for all A > 1. 



B. Finite Temperature 

Finally, let us present results for the T-dependence D{T) as 
evaluated using the relation (|6]l, again restricting our analysis 
to zero magnetization s = and systems with odd L (Fig. 7). 
It should be realized that numerical results at low T < 0.5 are 
more susceptible to finite-size effects since very small number 
of MB levels effectively participate in D{T) and the crucial 
contribution comes from the ground state eo(</i). Still, in spite 
of some discrepancies at low T < 0.4 the overall agreement 
with the TBA result lll3n is reasonable. Another conclusion is 
that the extended high-T behavior, i.e. D — C/T is followed 
very accurately down to quite low T > 0.5 in the whole range 
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Figure 7. (Color online) Finite-T Drude weight D(T) for A < 1 
at magnetization s = as calculated numerically using ED with 
1/ = 15 — 21 and finite-size scaling (full line with dots), extrapo- 
lating the high-T numerical result, i.e. D = C /T (thin lines), and 
Thermodynamic Bethe Ansatz result (dots) from Ref. Ll3.1 . 



A < 1. While the ground state value Dq is quite reliable in 
the intermediate window < T < 0.5 results are sensitive to 
finite-size effects so we cannot give a firm conclusion on pos- 
sible nonanalytical low-T behavior as predicted in Ref. |131 . 



IV. CONCLUSIONS 

In conclusion, we have shown that numerical evaluation of 
the Drude weight (spin stiffness) D{T) within the anisotropic 
Heisenberg model can lead to more controlled and converged 
results if performed in a canonical ensemble, at fixed Sz 
(number of particles N). Breaking of the particle-hole sym- 
metry by using systems with odd L is also helpful and is ad- 
vantageous over usually studied systems with even L. Our 
study is mostly concentrated on the high-T limit which should 
be anyhow quite generic for the whole regime T > 0. Results 
obtained at zero magnetization s = using the finite-size scal- 
ing confirm the change of character of D{T) at A = 1, i.e., 
they are compatible with the i:'(r) = for A > 1. While 
at s = within the majority of the regime A < 1 there are 
no evident problems with the scaling 1/L of D{T) we have 
traced the irregularities at A < 0.5 back to the emergence 
of finite-size low-w contribution in areg {'^) which can lead to 
a finite correction 5C in the thermodynamic limit L — > oo. 
Taken the latter into account, we find a very good agreement 



with the TBA result 1 13], in this way possibly eliminating (or 
at least restricting) some recently expressed questions regard- 
ing its validity. 

High-T normalized Drude weight D* away from s = 
shows a systematic and smooth variation with s towards the 
limiting values D* = 1 for s = ±1/2 as well as in XY limit 
A = 0. In the Ising regime A > 1 (in particular for large 
A > 2) the variation C = TD{s) is very well reproduced 
with the Mazur inequality overlap with the conserved energy 
current je, in very contrast to the XY-type regime A < 1. 

Results for the T-variation D{T) reveals that even quanti- 
tatively the high-T result D — C/T remains valid in a wide 
regime, i.e., generally for T > 0.5. While small-system re- 
sults allow also for a reliable scaling for Dq — D{T = 0) 
at s = 0, the finite-size effects are rather hard to avoid in the 
window < T < 0.5 and other methods beyond the ED are 
needed to investigate in more detail this regime. 
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